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CONVERTED-WAVE PRESTACK IMAGING AND 
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Summary 

A new technique for converted-wave (C-wave) prestack imaging and velocity analysis is 
presented. This technique, called pseudo-offset migration (POM) ? differs from 
equivalent-offset migration (EOM) (Bancroft et al., 1998) in how the migration mapping 
offset is defined. In this paper, we show that POM is less dependent on the initial 
velocity than EOM and can produce better images and C-wave velocities than P-S DMO 
plus poststack migration. 

Introduction 

Converted-wave processing is more difficult than pure-mode processing. A simple 
common-conversion-point stack requires an input trace to be grouped to different 
conversion-point positions. Velocity analysis becomes problematic even for a plane- 
layered medium. Migration is so critically dependent on both P-wave and S-wave 
velocities that a robust, efficient, and easy-to-use prestack migration method is needed 
for the C-wave processing. 

Equivalent-offset migration (EOM) is a prestack time migration method proposed by 
Bancroft et al. (1998). EOM performs prestack migration in two steps. The first step is 
the migration mapping. It maps each input sample to a common-scatter-point and 
accumulates it at an equivalent-offset (Figure 1). The equivalent-offset (h e ) is defined as 
a distance between the common-scatter-point (I) and a location of a collocated model- 
source (S e ) and model-receiver (Re) so the travel time from the model-source to the 
model-receiver equals to the input travel time from the source (S) to the receiver (R). 
The second step is the migration stack, which completes the migration by summing the 
common-scatter-point gather or EOM gather along the normal-mo veout travel times. 
EOM is claimed to be computationally efficient and weakly velocity-dependent. 
However, in C-wave processing, P-S EOM (Wang et al, 1996) depends on the initial 
velocity even for a single constant-velocity flat layer. A new scheme of migration 
mapping is proposed in this paper. We show that this new technique, called pseudo- 
offset migration (POM), is less velocity-dependent than EOM and produces better 
images than P-S DMO plus poststack migration (Harrison, 1 992). 

Pseudo-offset migration 

A pseudo-source (S p ) and a pseudo-receiver (R p ) are defined in Figure 2 such that the 
following two conditions are satisfied. 1) The total travel time from S p to D and from D 
to R p equals to the total travel time from S to D and from D to R. 2) The ray from S p to 
D has a common ray parameter p as the ray from D to R p . The distance x between S p and 
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R p is called the pseudo-offset. Unlike EOM, POM maps the input geometry to a new 
geometry defined by the pseudo-offset After the migration mapping, the migration stack 
or the C-wave velocity analysis can be performed on POM gathers based on the travel 
time from S p to R p . The travel time is exactly the same as the C-wave reflection travel 
time for the plane-layered media, which can be approximated by a high-order equation 
controlled by the C-wave RMS velocity (Thomsen, 1999). POM further distinguishes 
from EOM in that the pseudo-source and the pseudo-receiver can not be collocated 
together for C-waves because for the collocated source and receiver the C-wave ray 
parameters (p ep and p es in Figure 1) are not constant. C-wave velocity analysis using 
POM is independent of migration input velocities for the plane-layered media, but EOM 
is not. For dipping layers, POM mapping is dependent on the input velocities. However, 
C-wave velocity analysis becomes much easier after POM because the positions of 
dipping events are close to their true conversion-point positions. 

Figure 3 shows the velocity sensitivity test for both POM and EOM. All the gathers are 
moveout corrected with the correct model velocity. The model is a single layer with the 
constant C-wave velocity at 2121 m/s and the velocity ratio /q = 2. Figure 3a shows that 
if the input parameters are correct, both POM and EOM gather can be flattened with the 
model parameters. However, if the C-wave input velocity is wrong and even if y Q is still 
correct, EOM gather no longer has correct velocity information. The gathers in Figure 
3b are generated with an incorrect input velocity. Figure 3b shows that correcting the 
EOM gather using the model velocity can not flatten the gather. This indicates that the 
new velocity picked from the EOM gather depends on the input velocity. The POM 
gather shown in Figure 3b can be flattened with the model velocity. This indicates that 
the new velocity analysis of the POM gather does yield the correct velocity and it is 
independent of the input velocity. Even if both the input velocity and y$ may go wrong, 
figure 3 c shows the POM gather is less affected by the input parameters, but the EOM 
gather can no longer preserve the model velocity. 

Figure 4c shows a field data example of C-wave prestack time migration using POM. 
Comparing with P-S DMO plus poststack migration (Figure 4a), Figure 4c is well 
focused around the salt dome. Figure 4b shows a C-wave velocity spectrum after P-S 
DMO at the location A marked on Figure 4a. Figure 4d is the corresponding velocity 
spectrum after POM. Velocity analysis is much easier using Figure 4d than Figure 4b 
since the energy is well focused and is at the migrated position after POM. 

Conclusions 

POM is presented for C-wave prestack imaging. It is less sensitive to the velocity error 
than EOM. C-wave velocity analysis using POM is much easier than that using P-S 
DMO. POM is a useful technique for prestack time migration and velocity analysis. 
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Fig. 1. Geometry of equivalent-offset h e 
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Fig. 2. Geometry of pseudo-offset x 
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Fig. 3. Moveout corrected POM and EOM gathers. The model has the C-wave velocity 
at 2121 m/s and yo at 2. The input parameters used for migration mapping are 
(a) C-wave velocity 2121 m/s and yo 2 ? (b) C-wave velocity 1500 m/s and yo 2, 
and (c) C-wave velocity 1500 m/s andyo 3. 
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Fig. 4. Field data example: (a) P-S DMO followed by poststack migration, (b) C-wave DMO 

velocity spectrum, (c) prestack time migration using POM, and (d) C-wave POM velocity 
spectrum. 
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Converted-wave prestack time migration for isotropic and anisotropic media 

Weizhong Wang*, LongD. Pham, and Min Lou, PGS Seres Division, Houston, Texas, UX 



Summary 

We present a Kirchhoff prestack time migration technique for converted- waves (C-waves). This technique computes 
the C-wave diffraction travel time based on the exact travel time equations for v(z) VTI media. We demonstrate mat it 
produces more accurate image than migration algorithms using travel time approximations. For anisotropic media, this 
technique requires four parameters: C-wave interval velocity v c „ interval vertical-travel-time ratio y 0t , and Thomsen 
anisotropic parameters s x and For isotropic media, only v C[ and y 0l are required for the migration. We use pseudo- 
offset migration (POM) (Wang and Pham, 2001) to estimate C-wave velocities directly from the C-wave data. We 
show that C-wave velocity analysis becomes much easier on POM gathers because the positions of dipping events are 
closer to their true conversion-point positions and conflicting events are reduced after POM. 



Introduction 



Kirchhoff prestack time migration is often implemented with the hyperbolic travel time approximation, which is strictly 
correct for a homogeneous constant velocity medium. For multilayered v(z) media or vertical transversely isotropic 
(VTI) media, more sophisticated travel time equations are needed (Tsvankin and Thomsen, 1994; Li et al., 2001). The 
three-term approximation in VTI media introduced by Tsvankin and Thomsen (1994) provides sufficient accuracy for 
P-wave, but for SV-wave (abbreviated to S-wave), the area of validity of approximation is much limited. For C-wave 
reflection moveout, the three-term approximation can still be accurate up to the offset-to-depth ratio of 2.0 (Thomsen, 
1999; Yuan, 2000) . C-wave velocity analysis and common-conversion-point (CCP) binning using the three-term travel 
time approximation can produce accurate results since C-wave data is normally muted at large offset-to-depth ratios. 
However, for C-wave prestack migration, the travel time approximation is no longer acceptable for imaging high 
dipping events. The S-wave ray-path is longer than the P-wave ray-path when the scatter point is close to the source 
location. The migration radius-to-depth ratio could easily be larger than an acceptable level. In this paper, we compute 
C-wave diffraction travel time using exact travel time equations and show the improved dipping image by comparing 
with the migration using the three-term approximation of Tsvankin and Thomsen (1 994). 



Migration parameters that we use for computing the exact travel times are C-wave interval velocity va, interval 
vertical- travel-time ratio y 0l , and Thomsen (19S6) anisotropic parameters s l and The v Cf is computed from the C- 
wave short-spread moveout velocity V C2 using the Dix formula. The y 0i has to be determined by correlating 
corresponding P-wave and C-wave events from stacked sections or near-offsets of prestack gathers. This correlation 
process is also required by CCP binning and NMO correction (Thomsen, 1999). The only difference here is that we 
need to convert the vertical-travel-time ratio y$ to the interval ratio y 0l . The estimation of the anisotropic parameters 
(Alkhalifah and Tsvankin, 1995; Li and Yuan, 2001) is not addressed in this paper. Our main focus is how to estimate 
2^ accurate C-wave velocities. A pseudo-offset migration (POM) method (Wang and Pham, 2001) is presented in this 

UJ paper for this purpose. 



& Prestack Time Migration with Exact Travel Times 

O 

^ The diffraction travel time at a scatter point in Kirchhoff prestack time migration is the total travel time of P-wave (/ P ) 

down from the source to the scatter point and S-wave (t s ) up from the scatter point to the receiver: 

f = f P +f s . (1) 
For v(z) VTI media, we have the following equations: 

g p = ^ Ai ™< v PQi/ v ?i(P?> v ?Qi> v ^ £ i A ) ^ ( 2 ) 

i=l Pl V l(P?^P0n V S0t^ A) 

f s =^ ^SOi v S0iA&O?S> v P0i A ) ^ ^ 

1=1 ^-pMXPs>v ? < h >Vsoi>£iA') 

h ? = ^ P? v rt(P?> v ™> v sv> £ i A t >m At m ^ 
1=1 V*~ PpVwOjp.vpo^Vgo^e, ,8, ) 

^ = V Is^O^po^so^ A )yfcoiA*so t ^ ^ 
i=i V* ~ />s4(Ps> v PO;> v so;> f f A ) 



where h ? and h s are the horizontal distances from the source / receiver to the scatter point, p P and p$ are the P- and S- 
wave ray parameters, A/ POl and A?so, are the vertical-travel-times in the i-th layer, v POl and v S o t are the vertical velocities 
of the i-th layer, and v ?l and v S/ - are the group velocities. Based on the work of Thomsen (1986), we can derive the group 
velocities as a function of p? / v P0( , v SO i, s x and S t numerically. The ray parameters in equations (2)-(5) are not time- 
variant, and h ? and h s are constants for a given scatter point, so we can numerically solve equations (4) and (5) for p P 
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and p$, To compute travel times f P and t s , we need to know A2 P0[ , Az S0/ , v P0( -, v S0( , e ti and 8 V The A/ SOl - and v S0 ; values can 
not be measured or estimated directly from P- and C-wave data. We replace them by the C-wave vertical-travel-time 

A'co/^poi+A'soi* ( 6 ) 
and the interval vertical-travel-time ratio, that is the same as the vertical-velocity-ratio for the i-th layer, 

r 0 i = ^so, / A* pw = v po ( - / v so; ■ ( 7 ) 
With equations (6) and (7), At m can also be replaced by Afc 0 , and y 0i . The v POt parameter can be determined from P- 
wave short-spread moveout velocity v P2i , but the P-wave velocity is indexed at P-wave times and has to be mapped to 
C-wave times using y oi . The y 0l parameter is determined by correlating corresponding P-wave and C-wave events from 
stacked sections or near-offsets of prestack gathers. Thus, it will have some errors. P-wave velocity v P2l often has its 
own error. If the P-wave velocity error is propagated through y 0l to C-wave migration, it is much difficult to take it out 
without joint processing of both P- and C-wave data. It is better to use a C-wave velocity to control the C-wave 
migration. Using the definition of the C-wave interval velocity of Thomsen (1999), 

v 2 = jL + j4 (g) 

i+r 0l l+i/n, 

where 

v P 2 2i -v P 2 0l (l+2^), (9) 
Vs 2l =v s 2 0j (l + 2O, (10) 

we replace v P0/ and v SOi with v Cl) y Qu % and S v Thus the total diffraction travel time can be accurately computed with 
four parameters Vq, ftu 3 an ^ 4 at a given C-wave vertical-travel-time f co . 

In our equations, we do not have ^(Thomsen, 1999) which requires to know P-wave velocities. The v(z) effect in 
travel times is reflected in the interval velocity v c/ and the interval y Ql . The VTI effect is reflected in anisotropic 
parameters a, and S l directly. For isotropic v(z) media, equations (2)-(5) can be simplified and the diffraction travel time 
depends on only two parameters v Cf and y Q[ for a given t co . 

Figure 1 shows a synthetic example of C-wave prestack time migration for a single layer VTI model. The Kirchhoff 
migration with exact travel times computed using equations (l)-(5) is shown in Figure lb. Comparing with Figure lc 
which is the migration with travel time approximation of Tsvankin and Thomsen (1994), the migration with the exact 
travel times produces well-focused image for the large dipping section. 

The exact travel time computation does not increase the migration runtime much since for prestack time migration a 
small travel time table could be built at each scatter point for all relevant input traces. 



Pseudo-Offset Migration 



Bancroft et aL (1998) developed the equivalent-offset migration (EOM) method for Kirchhoff prestack summation. 
Wang and Pham (2001) modified it for converted-wave data by introducing a new offset, called pseudo-offset, and 
pointed out that migration with pseudo-offset (POM) is less dependent of the initial velocity than EOM. We use POM 
here to estimate the C-wave velocity for v(z) VTI media. 

POM performs the Kirchhoff summation in two steps. The first step is the migration mapping. For each given scatter 
point D, POM maps each input sample from the survey geometry, defined by the source S and the receiver R in Figure 
2, to a new geometry, defined by a pseudo-source S* and a pseudo-receiver R*. There are two conditions the new 
geometry must be satisfied. 1) The total travel time from S* to D and from D to R* must be equal to the total travel 
time from S to D and from D to R. 2) The ray from S* to D and from D to R* must have a common ray parameter p. 
With these two conditions, we can compute the ray parameter p by solving the following equation 

h (p) + 's (P) = h (P? )+*s(P&)> ( 12 ) 
and then compute the pseudo-offset, that is defined as 

* = ApO>)+AsO>). ( 13 > 
The input samples mapped to the same pseudo-offset are accumulated without any time shift so a common-scatter point 
gather is built with different pseudo-offsets, which is called the POM gather. The second step in POM is the migration 
stack, which completes the migration by summing the POM gather along the moveout travel times, 

t = t P (p) + t s (p), (14) 
where p is computed from equation (13) for each given pseudo-offset x. Other parameters used for computing the 
moveout travel times in equation (14) are v C ;, % w £i and 6 V Assuming y 0i , s x and S t without errors, the reflection moveout 
of the plane layers is not altered by the migration mapping even if the initial vc* may be wrong. Thus, POM gathers can 
be used to estimate correct v Cf for plane-layered media. Figure 3a shows an example of a POM gather from a flat 
section in the VTI model of Figure la. The POM gather is created using a wrong initial velocity at 2000 m/s and is then 
corrected with correct model parameters. The flattened event in Figure 3a indicates that we can estimate correct model 
parameters from it even if the initial C-wave velocity is off by almost 30%. However, the direct inversion of C-wave 
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interval velocity v Cl using the moveout equation (14) is difficult. One can use the following approximation to scan for 
the C-wave short-spread moveout velocity V C2 (Thomsen, 1999; Li and Yuan, 2001): 



-), A- 



(15) 



and then compute vc using Dix formula. Figure 3b shows that the moveout correction using equation (15) is acceptable 
only in a limited offset range. Equation (15) is an approximation for a single-layered isotropic medium. To improve its 
accuracy for velocity analysis of v(z) VTI media, we suggest to compute the residual moveout between equation (14) 
and (15) using the initial parameters, correct POM gathers with the residuals, and then scan POM gathers with equation 
(15) for correct velocity Kc2- Figure 3c shows that this procedure can produce almost an identical result as the exact 
travel time equation (14). 

Figure 4c shows a field data example of isotropic C-wave prestack time migration using POM. Comparing with C- 
wave DMO plus poststack migration (Harrison, 1992) (Figure 4a), Figure 4c is well focused around the salt dome. 
Figure 4b shows a C-wave velocity spectrum after C-wave DMO at the location A marked on Figure 4a. Figure 4d is 
the corresponding velocity spectrum after POM. Velocity analysis is much easier using Figure 4d than using Figure 4b 
since the energy is well focused and is at the migrated position after POM. 



Conclusions 



Travel time computation must be accurate enough for imaging high dipping events of converted-wave data. We derive 
the exact travel time equations for v(z) VTI media using four parameters: v c „ fou fy, and 5 V We present pseudo-offset 
migration to perform Kirchhoff prestack summation and C-wave velocity analysis. We demonstrate that even if the 
initial velocity is wrong, POM gathers still preserve correct moveout for estimation of C-wave velocities of flat layers. 
For complex structures, C-wave velocity analysis on POM gathers is much easier than on CCP or DMO gathers. 
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Fig. 1. Converted- wave prestack time migration of a VTI 
model, (a) VTI model of Taylor sandstone, (b) 
Migration with exact travel times, (c) Migration with 
travel time approximation. 
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Fig. 2. Geometry of pseudo-offset x 
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Fig. 3. Moveout corrected POM gathers with the correct C-wave 
velocity (2830 m/s) for the Taylor sandstone model. The 
initial velocity used in migration is wrong at 2000 m/s. The 
travel time equations used for moveout correction are (a) 
the exact moveout equation (14), (b) approximation 
equation (15), and (c) equation (15) but with residual 
moveout correction using the initial velocity. 



Fig. 4. Field data example: (a) DMO followed by poststack 

migration, (b) C-wave DMO spectrum, (c) prestack time 
migration using POM, and (d) C-wave POM velocity 
spectrum. 



ATTOmmY DOCKET NO. 30290US 



